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Abstract. In this paper we investigate the asymptotics of forward-start options and the forward implied 
volatility smile in the Heston model as the maturity approaches zero. We prove that the forward smile for out-of- 
the-money options explodes and compute a closed-form high-order expansion detailing the rate of the explosion. 
Furthermore the result shows that the square-root behaviour of the variance process induces a singularity 
such that for certain parameter configurations one cannot obtain high-order out-of-thc-moncy forward smile 
asymptotics. In the at-the-money case a separate model-independent analysis shows that the small-maturity 
limit is well defined for any Ito diffusion. The proofs rely on the theory of sharp large deviations (and refinements) 
and incidentally we provide an example of degenerate large deviations behaviour. 
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p i . 1. Introduction 

Consider an asset price process (e Xt ) 4>0 with = 0, paying no dividend, defined on a complete filtered 
'-4_ — 
I 1 probability space (f2, J 7 , (J-'t)t>o, P) with a given risk-neutral measure P, and assume that interest rates are 

zero. In the Black-Scholes-Mcrton (BSM) model, the dynamics of the logarithm of the asset price are given by 

■ dX t = — A(7 2 di + adW t , where a > represents the instantaneous volatility and W is a standard Brownian 

t> i motion. The no-arbitrage price of the call option at time zero is then given by the famous BSM formula [TUll4Tj : 

\Q ■ Cbs(t, k, a) := E (c Xt — e fe ) + = N {d + ) — c k N (cL), with d± := --^p±\o^pr ', where N is the standard normal 

distribution function. For a given market price C obs (r, k) of the option at strike e fc and maturity r, the spot 

, implied volatility a T (k) is the unique solution to the equation C obs (r, k) — Cbs(t, k, a T (k)). 

, For any t, r > and k £ R, we define as in [9l|40] a forward-start option with forward-start date t, maturity r 

and strike e k as a European option with payoff (e x ^ ' — e k \ where we define Xt^ := Xt+ T — Xt pathwise. By 



the stationary increment property, its value is simply Cbs( t , k,a) in the BSM model. For a given market price 
"^^"j , C° hs (t, r, k) of the option at strike e fc , forward-start date t and maturity r, the forward implied volatility smile 
d ' °t,T(fc) is then defined (see also [HI [20]) as the unique solution to C obs (<, r, k) = C-qs( t i k, a t T (k)). The forward 
smile is a generalisation of the spot implied volatility smile, and the two are equal when t = 0. 

Asymptotics of the spot implied volatility surface have received a large amount of attention over the past 
decade. These results have helped shape calibration methodologies based on arbitrage-free approximation of 
the spot smile in a large variety of models. Small-maturity asymptotics have been studied by Berestycki- 
Busca-Florent [8] using PDE methods for continuous time diffusions and by Henry-Labordere [28] using heat 
kernel expansions. Forde et al. |18) and Jacquier et al. |32| derived small and large-maturity asymptotics in 
the Heston model and affinc stochastic volatility models with jumps using large deviations and saddlcpoint 
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methods. Wing asymptotics (as \k\ /" oo) have been studied by Lee [39] and further extended by Benaim 
and Friz [SJ [B] and in [T31 HU [57] . Fouque et al. [ID] have used perturbation techniques in order to study slow 
and fast mean-reverting stochastic volatility models. Small-noise expansions have been studied by Osajima [33j 
and Takahashi J35] using Watanabe expansions and in Deuschel et al. [T3] using Laplace method on Wiener 
space. See [26l [48] for a general overview of implied volatility asymptotics in stochastic volatility models. In 
exponential Levy models implied volatilities of out-of-the-money options explode as the maturity tends to zero, 
while the implied volatility for at-the-money options converges to the volatility of the diffusion component as 
the maturity tends to zero. Small-maturity asymptotics for models with jumps (including Levy processes) have 
been investigated in [2] [3] |46 ] |43 ] |42 ] [TB ] . 

However, these asymptotics do not provide any information on the forward smile or forward-start-based 
payoffs, such as cliquets and forward-start options and the literature on this topic is sparse. Glasserman and 
Wu [23] introduced different notions of forward volatilities for forecasting purposes, Keller-Ressel [33] studies 
a very specific type of asymptotic (when the forward-start date t becomes large), and empirical results have 
been carried out in [pj 1231 lllj . In |15j the authors empirically studied the forward smile in Sato models and 
ran comparisons with a suite of models including Heston and local volatility models for forward smile sensitive 
products such as cliquets. Recently, small and large-maturity forward smile asymptotics were derived in [3T] 
for a general class of models including the Heston model. However, the results in [31] only apply to the so- 
called diagonal small- maturity regime, i.e. the behaviour (as e tends to zero) of the process (X^) £ >o- The 
conjecture, stated there, is that for fixed t > the Heston forward smile (corresponding to X^) explodes to 
infinity (except at-the-money) as r tends to zero. 

In this paper we confirm this conjecture and give a high-order expansion for the forward smile. The main 
result (Theorem 14. ip is that the small-maturity Heston forward smile explodes according to the following 
asymptotic: of T (fc) = v (k,t)r~ 1 / 2 + vi(k,t)r~ 1 / 4 + o (t -1 / 4 ) for k € R* and t > as r tends to zero. Here 
vo(-,t) and v\(-,t) are even continuous functions (over R) with vo(Q,t) = i>i(0,i) = and independent of the 
Heston correlation. In the at-the-money case {k = 0) a separate model-independent analysis (Lemma 14.31 and 
Theorem 14. 4[) shows that the small-maturity limit is well defined and lim T ^o 0t. T (O) = ^WVt) holds for any 
well-behaved diffusion where Vt is the instantaneous variance at time t. This exploding nature is consistent 
with empirical observations in |llj and the diagonal small-maturity asymptotic from [31] . 

The paper is structured as follows. In Section [5] we introduce the notion of a forward time-scale and charac- 
terise it in the Heston model. In Section [3] we state the main result on small- maturity asymptotics of forward- 
start options in the Heston model. Section @] tackles the forward implied volatility asymptotics: Section 14.11 
translates the results of Section [3] into out-of-the-money forward smile asymptotics, and Section B~2l presents a 
model-independent result for the at-the-money forward implied volatility. Section[5jprovides numerical evidence 
supporting the asymptotics derived in the paper and the main proofs are gathered in Section [6] 

Notations: E shall always denote expectation under a risk-neutral measure P given a priori. We shall refer 
to the standard (as opposed to the forward) implied volatility as the spot smile and denote it a T . The forward 
implied volatility will be denoted at. T as above. In this paper we will always assume that forward-start date is 
greater than zero (t > 0) unless otherwise stated. 
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2. Forward time-scales 

In this section we introduce the notion of a forward time-scale and characterise it in the Heston model. In 
the BSM model the time-scale is given by h(t) = t, which is related to the quadratic variation of the driving 
Brownian motion. For diffusions (such as Heston) the (spot) time-scale is the same, which implies that the spot 
smile has a finite (non-zero) small-maturity limit. In the forward case, this however no longer remains true. 
Stochastic volatility models (eg. Heston) exhibit different time-scales to the BSM model leading to different 
asymptotic regimes for the forward smile relative to the spot smile. As we will show below (Lemma 12. 3[) , all 
re-scalings of the Heston model lead to limiting mgf 's that are all zero on their domains of definition. But the 
forward time-scale is the only choice that leads to the limiting mgf being zero on a bounded domain. This is 
one of the key properties that allows us to derive sharp large deviation results even though at first sight this 
zero limit appears trivial and non-consequential. We define the forward moment generating function by 

(2.1) A^(w,a) := alogE U"*™ 'A , for all u £ V t>T , 

where T> t>T := {ti £ R : lA^^a)! < oo}. With this definition the domain 2? tjT will depend on a, but it will be 
clear from the context which choice of a we are using. 

Definition 2.1. We define a (small-maturity) forward time-scale as a continuous function h : K+ — > R+ such 
that lim^o h(r) = and A(it) := lim^o At (u, h(r)) produces a non-trivial pointwise limit. We shall say that 
a (pointwise) limit is trivial if it is null on M or null at the origin and infinite on R* . 

Remark 2.2. 

(i) In the BSM model the forward time-scale is h(r) = r. 

(ii) A forward time-scale may not exist for a model. For example, consider exponential Levy models with 
bounded domain for the Levy exponent. The only non-trivial limit occurs when h = 1, which does not 
satisfy Definition 12.11 and so a forward-time scale does not exist. 

(iii) If (Xr^) T >o satisfies a large deviations principle [T^l Section 1.2] with speed h and assuming further some 
tail condition (see Theorem 4.3.1]), then h is the forward time-scale for the model by Varadhan's 
lemma. 

(iv) Diffusion models have the same spot time-scale {t = 0) as the BSM model, namely h(r) = r (see for 
example [5]). This is not necessarily true in the forward case as we will shortly see. 

In the Heston model, the (log) stock price process is the unique strong solution to the following SDEs: 

dX t =-^V t dt+^V t dWt, A o = 0, 

( 2 - 2 ) dV t = n(8-V t )dt + ZVV t dZ t , V o = v>0, 

d(W,Z) t =pdt, 

with k>0, £ > 0, 8 > and \p\ < 1. The Feller SDE for the variance process has a unique strong solution by 
the Yamada-Watanabc conditions [331 Proposition 2.13, page 291]). The X process is a stochastic integral of the 
V process and is therefore well defined. The Feller condition, 2k8 > £ 2 , ensures that the origin is unattainable. 
Otherwise the origin is regular (hence attainable) and strongly reflecting (see [3J] Chapter 15]). We do not 
require the Feller condition in our analysis since we work with the forward mgf of X which is always well defined. 
In order to characterise the Heston forward time-scale we require the following lemma, proved in Section 16.11 
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Lemma 2.3. Let h : R + — > R + be a continuous function such that lim^o M T ) — and a £ R5_. The following 
limits hold for the Heston forward mgf as t tends to zero with (3 t defined in (|6.3p : 

(i) If h(r) = a^fr then \hn T \^o A^\u, h(r)) = 0, for all \u\ < aj\fB~t o-nd is infinite otherwise; 

(ii) if y/r/h^r) oo then lim^o At \u, h(r)) = 0, for u = and is infinite otherwise; 
(Hi) ify/r/h(r) \ then linv^o A T t} (u, h(r)) = 0, for all u £ R. 

As it turns out all limits are zero on their domains of definition, but using h(r) = y/r produces the only 
(up to a constant multiplicative factor) non-trivial zero limit. It follows that t i— > yfr is the Heston forward 
time-scale. Let now A : T>\ = (—l/y/Jfl, l/^ffit) — > R be the pointwise limit (/3t in (|6.3|) ) from Lemma \2.i\ i.e. 
satisfying A(u) = for u £ Dk and infinity otherwise. Further we define the function A* : R — > R + as the 
Fenchel-Legendre transform of A: 



(2.3) 



A*(fc) := sup {uk — A(it)} , for all k £ 
uev A 



Lemma 2.4. The function A* defined in (|2.3|) is characterised explicitly as A*(fc) = \k\/y/]3t for all k £ R. 
Proof. Clearly A*(0) = 0. Now suppose that k > 0. Then A*(fc) = sup ugX , A {uk} = kj \fB~ t . A similar result 



holds for k < and the result follows. 

3. Small-maturity forward-start option asymptotics 



□ 



In this section we state the main result on small-maturity forward-start option asymptotics. First we need to 
define a number of functions. All functions below are real- valued and defined on R*. We recall that sgn(u) = I 
if u > and -1 otherwise and fit is defined in (|6.3|) . 



(3-1) 

where 
(3.2) 

(3.3) 
(3.4) 



(3.5) 



sgn(fc) 
a (k) := — oi(fc) 



a (fc)v^ Kt/2 



a 3 (k) 



2P t a\{k) 



k\f3, 



1/4 



K e Si(oo(fc)) 



^v 2 



ev(3 t c Kt (\k\ePl Bi(ao(fc)) - fc£ 2 .§iMfc)) - «*) + (2K0(3 t c Kt ) 2 



B^u) :=^[u 2 pt-2 



a(k) :-- 



2^(T Kt l 2 
e (k) 3 / 2 ' 



•(*) 



a 2 (k) 



k6 



\k\ev& 7 



e (k) := -2ai(fc)/a (fc), 
ei(fc) := -2f3 t r(k), 

e 2 (k) := -2/3 f (ai(fc)a a (fc) + a (fc)a 3 (fc) + Oi(fc)^(ao(fc)) 



^i(A) 



ao(k)ve 



eg(fc) + ao(fc)A [3ai(fc)e (fc) - 2a (fc)ei(fc)] 



4a (fc)^ t e- K 7e4(fc) 

(4flo(ft)A[3ao(fc)ei(A) - 4 ai (fc)e (fc)] - 5e 2 (fc)) , 



we 

8v/3 t e- Kt /e 5 (k), 



vc 



2el{k) 



j(k) - e (fc)e 2 (fc) 
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2a (fc)ai(fc)e (fc)ei(fc) + 2e 2 (k)r(k) , 
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cj>l(k) := ip 3 {k) - t(jo{k)ipi{k) - 2 m — , 

4> c 2 (k) :=-V?(fc)/2, 



^ (u\ , 2^ ei (fc) <^(fc) 3^(fc) , 1501 (fc) 

(3.7) Zl (k) := Mk) - a 3 (k)k - ^^), Px(*) := eo(*) + ^ + ^V" + 

/ c (fc) :=2|oi(*)fc|, c 1 (fc):=^fa (fc)ai(fc)--|i% > )-a 2 (A ; )fc ) 

(3.8) < e (fc) V 2/3 t e (k)J 

{ cb(*;) :=e (fc)- 2Ke/ « 2 , c 3 (fc) :=Z!(fc)+ Pl (fc). 

We now state the main result of the section, i.e. an asymptotic expansion formula for forward-start option 
prices as the remaining maturity tends to zero. The proof is given in Section [6.41 



Theorem 3.1. The following expansion holds for forward- start option prices for all fcel* as r tends to zero: 
E 



c A - - e fc 



= (l-c fc )l {fe <o } 

/ A*(k) c (fc) /IN ,\ /3 t r( 7/8 - eK/(2 « 2) )c 2 (fc) / ,,, 1M / 1/4 

where A* is characterised in Lemma \2.4\ CQ,...,ca in (|3. 8|) . <r in (|3.3p and /3 t m (|6.3[) . 
Remark 3.2. 

(i) We have A*(k) > and co(fc) > for all fceM*. Also note that A* is linear as opposed to being strictly 
convex in the BSM model, see Lemma I3.4I below. 

(ii) The forward time-scale y/r results in out-of-the-money forward-start options decaying as r tends to zero 
at leading order with a rate of exp (— 1/y/r) as opposed to a rate of exp (— 1/r) in the BSM model. 

(iii) The fact that the limiting forward mgf is non-steep (trivially zero on a bounded interval) results in a 
different asymptotic regime for higher order terms compared to the BSM model. In particular we have a 
t 1 / 4 dependence as opposed to a r dependence in the BSM model and the introduction of the parameter 
dependent term r < - 7 / 8_eK /' 2 ^ ^. The implications of this parameter dependent term for forward-smile 
asymptotics will be discussed further in Remark I4.2f iv) . 

(iv) The asymptotic expansion is given in closed-form and can in principle be extended to arbitrary order using 
the methods given in the proof. 



As an immediate consequence of Theorem 13.11 we have the following corollary, which provides an example of 
a family of random variables for which the limiting moment generating function is zero (on its effective domain) 
but a large deviation principle still holds. This is to be compared to the Gartner-Ellis theorem [TU Theorem 
2.3.6] which requires the limiting mgf to be at least steep at the boundaries of its effective domain for an LDP 
to hold. 

Corollary 3.3. ( Xf ) satisfies an LDP with speed \fr and good rate function A* as r tends to zero. 



Proof. The proof of Theorem 13 . 1 1 holds with only minor modifications for digital options, which are equivalent to 
probabilities of the form P (x^ < k) or P (x { T t] > fc) . One can then show that lim rN0 ^/r\og¥ < kj = 
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— inf{A*(x), x < k}. Note that of course this infimum is null whenever k > 0. Consider now an open interval 
of the real line of the form (a,b). Since (a, b) = (—00,6) \ (—00, a], then by continuity of the function A* and 
its properties given in Lemma 12.41 we immediately obtain that 



lim Jt 

t\0 



log P fx® e (a,b)) = - inf A*(x). 

V / x£(a,b) 



Since any Borel set of the real line can be written as a (countable) union / intersection of open intervals, the 
corollary follows from the definition of the large deviations principle, see [T2j Section 1.2]. □ 

In order to translate the forward-start option results into forward smile asymptotics we require a similar 
expansion for the BSM model. The following lemma is a direct consequence of [3TJ Corollary 2.9] and the proof 
is therefore omitted. 

Lemma 3.4. In the BSM model the following expansion holds for all k £ R* as r tends to zero: 

,k/2-k 2 /(2a 2 T) 



E 



(e^'-e*)' 



l-e k )l {k<0} + 



k 2 V2^ 



4. Small-maturity forward smile asymptotics 



4.1. Out-of-the-money forward implied volatility. We now translate the small-maturity forward-start 
option asymptotics into forward smile asymptotics. Define the functions Vi : R* x M*, — > R (i = 0, 1, 2, 3) by 

fc 2 VJTt\k\ 2cg(k)v*(k,t) e^«/ 2 A 1/4 74fe| 
Mk,t) ^ — = -^-, Mk,t):= - 2 = , 

(41) ll(k A Hik.t) ( e^c 2 {k)[3 t e \ v*(k,t) v\{k,t) 

v 3 (k,t) := !^( 2c3 (fcMM) - 3 Vl (k,t)) + (2v 2 (k,t) - ^M) , 

with A* characterised in Lemma [2.41 Co, . . . , C3 in (|3.8p . a in (|3.3|) and /3 t in (|6.3|) . On R*, A*(fc) > and so 
vo(k,t) > 0. Further co(fc) > and so vi(k,t) > 0. Also C2(k) > and a(k) > so that v 2 is a well defined 
real- valued function. The following theorem — proved in Section 16.41 — is the main result of the section. 



Theorem 4.1. The following expansion holds for the forward smile for all k € R* as r tends to zero: 

2 



ff 2 = I 7-1/2 7-1/4 ^ T l/4 



Remark 4.2. 

(i) Note that Vo(k,t) and v±(k,t) are strictly positive for all fc e 1*, so that the Heston forward smile blows 
up to infinity (except ATM) as r tends to zero. 

(ii) Both vo(-,t) and vi(-,t) are even functions and correlation-independent quantities so that for small matu- 
rities the Heston forward smile becomes symmetric (in log-strikes) around the ATM point. Consequently, 
if one believes that the small-maturity forward smile should be downward sloping (similar to the spot 
smile) then the Heston model should not be chosen. This small-maturity 'U-shaped' effect for the Heston 
forward smile has been mentioned qualitatively by practitioners; see [11) . 
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(iii) We use the notation / ~ g to mean f jg — 1 as r — > 0. Then in Heston we have a\ T ~ v / S"l^l/(2v / ^ : ) an d 
in exponential Levy models with Levy measure v satisfying supp^ =lwe have of r ~ — /c 2 /(2rlogr) 
Page 21]. Wc therefore see that the small-maturity exponential Levy smile blows up at a much quicker 
rate than the Heston forward smile. 

(iv) We have limfc_>.o «o(fc, t) = vo(0,t) = and lim; c _j.o v\(k 7 t) = vi(0,t) = 0. Higher-order terms are not 
necessarily continuous at k = 0. For example (when Ak6 = £ 2 ) we have lim/ c _ i .o V2(k,t) = +oo. 

(v) The at-the-money forward implied volatility (fc = 0) asymptotic is not covered by Theorem 14.11 and a 
separate analysis is needed for this case (see Section |4"T2"|) . In particular the proof fails since in this case the 
key function m*(0) (defined through equation (|6.13p ) does not converge to a boundary point, but rather 
to zero as r tends to zero (see the proof of Lemma l£>.3[) . 

(vi) It does not make sense to consider the limit of our asymptotic result for fixed k G M* as t tends to zero 
since for t = using the forward time-scale h(r) = yfr will produce a trivial limiting mgf and hence none 
of the results will carry over. The time scale in the spot case is h(r) = r; see |17j . Our result is only valid 
in the forward (not spot) smile case. 

(vii) As seen in the proof, due to the term r 7 / 8 ~ eK /' 2 ^ ) in the forward-start option asymptotics in Theorem l3.ll 
one can only specify the small-maturity forward smile to arbitrary order if 4k8 = £ 2 . If this is not the case 
then we can only specify the forward smile to order ©(I/t 1 / 4 ). Now the dynamics of the Heston volatility 
o t := \fVt is given by da t = ( 4K ®~f - dt + |dW t , with a = y/v. If we set 4k0 = £ 2 the volatility 
becomes Gaussian, which corresponds to a specific case of the Schobel-Zhu stochastic volatility model [45] . 
So as the Heston volatility dynamics deviate from Gaussian volatility dynamics a certain degeneracy occurs 
such that one cannot specify high order forward smile asymptotics in the small-maturity case. Interestingly, 
a similar degenerary occurs when studying the tail probability of the stock price. As proved in |13j . the 
square- root behaviour of the variance process induces some singularity and hence a fundamentally different 
behaviour when Ak9 £ 2 . 

4.2. At-the-money forward implied volatility. The analysis above excluded the at-the-money case k = 0. 
Wc show below that this case has a very different behaviour and can be studied with a much simpler machinery. 
In this section, we shall denote the future implied volatility at(k,r) as the implied volatility corresponding to 
a European call/put option with strike e fc , maturity r, observed at time t. We first start with the following 
model- independent lemma, bridging the gap between the at-the-money future implied volatility cr t (0, r) and the 
forward implied volatility <7t. T (0). Note that a similar result — albeit less general — was derived in [3"Tj . 

Lemma 4.3. Let t > 0. Assume that there exists n € N* such that the expansion <7t(0, r) = Y^j=o tT i(^) rJ +° ( T ™) 
holds and that Eq (ct j(t)) < oo for j = 0, ...,n. Then the at-the-money forward implied volatility satisfies 

^,r(0)=E?= O E0 (<X^)H + 0(T«). 

Proof. In the Black-Scholcs model, we know that for any t > 0, r > 0, the price at time t of a (re-normalised) 
European call option with maturity t + t is worth C ( BS (k,r,a) := E (St+ T / St — e fc ) \Ft , and its at-the-money 
expansion as the maturity r tends to zero reads (see |18l Corollary 3.5]) 



Cf(0,T,o-) 



(a 3 / 2 )J. 



I^K V 24 

For a given model, we shall denote by Ct(fc, r) the price at time t of a European call option with pay- 
off (St+ T /St — e k ) at time t + r. The future implied volatility at{k,r) is then the unique solution to 
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Cf s (k, r, <xj(fc, r)) = Ct(k, t). At the money k = 0, we obtain the following expansion for short maturity t: 
(4-2) Q(0,r) = -i= ^ t (0,r)^- +o(a t 3 (0,r)r 3 / 2 )) . 

Now, the forward-start European call option (at inception) in the Black-Scholcs model reads 

St + T 



C^{k,r,a) :=E 



5, 



t+T 



5 t 



E < E 



1^ 



■■N~(d+{(T, t)) - e k M{d-{o-, r)), 



where d±(er, r) := (— k ± <t 2 t/2)/ (cry/r). For a given model, let us denote by C(k,t,r) the price of a forward- 
start European call option. By definition of the forward implied volatility c t]T (fc), we have C{k,t,r) = 
C^ s (k,T,a t . T (k)). For at-the -money k = 0, it follows that 



C(0,t,r) = C BS (0,r,at, T (0)) =^(d+K T (0),r)) - AT(d_(a t , r (0),r)). 



Using the expansion 



(4.3) 
Recall now that 
(4.4) C(k,t,r) :=E 



we obtain 



C(0,t,r) = 



q t , T (0) 3 r 3 / 2 
0t,r(O)VV + ° 



.3/2 



t+T _ c fc 



St 



= En 



St + T 



E (C t (fc,T)). 



Under the assumptions of the lemma, with n = 1, suppose as an ansatz that the corresponding forward smile 
satisfies ot, T (0) = X^=o^j W 1 ""' + °( r )- Combining this with (|4.2j) and (|4.3[) . we find that <fj(<) = E (<Jj(t)) for 
j = 0, 1. The higher-order terms for the expansion can be proved analogously and the lemma follows. □ 

We now apply this to the Heston model. Recall the definition of the Kummer (confluent hypergeometric) 
function M : C 3 -> K: 

(«)» z n 



M(a,ft,z) := J2 



n>0 



/^0,-l,... 



where the Pochhammcr symbol is defined by (a) T 
p > -2k0/£ 2 and t > we define 



a (a + 1) ■ • • (a + n — 1) for n > 1 and (a)o = 1- For any 



(4.5) 



A(t,p) := 2 p $ , exp 



2/8* 



_v 12 11 m I 1- n. 



vc 



p, 



r(2 K e/e) 

with Pt defined in (|6.3p . This function is related to the moments of the Feller diffusion (see [131 Theorem 
2.4]): for any t > 0, E [Vf] = A(t,p) if p > —2n8/£ 2 and is infinite otherwise. Note in particular that 
lim t N^ A(t,p) = v p (see [TJ 13.1.4 page 504]). The Heston forward at-the-money volatility asymptotic is given 
in the following theorem. 

Theorem 4.4. The following expansion holds for the forward at-the-money volatility as r tends to zero: 

ifine<e, 



A(t,l/2) + o(l), 
^ ( ° ) = < A(M/2) + A( *'- 1/2) 



24 VP V 8 



(p£- 2k)t + o(t), if4 K 0>^ 2 



Remark 4.5. 

(i) As opposed to the out-of-the-money case, the small-maturity limit here is well defined. 

(ii) Combining Lemma T4.3I and [5], lim T ^o vt,r (0) = E(y / Vt") holds for any well-behaved stochastic volatility 
model (S,V). 
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(iii) The proof does not allow one to conclude any information about higher order terms in Heston for the case 
4k9 < £ 2 . A different method would need to be used to compute higher order asymptotics in this case. 

Proof of Theorem \4-4\ In Heston we recall ( [181 Corollary 4.3] or [22 Corollary 3.3]) the asymptotic of (0,r) = 
V t + (if + (P 2 ~ 4 ) + - 2k )) t + o(t), and so for small r we have (T t (0,r) = a (t) +ai(t)r + o(r), with 



We first compare the true Heston forward smile and the asymptotics developed in the paper. We first 
calculate forward-start option prices using the inverse Fourier transform representation in |38[ Theorem 5.1] 
and a global adaptive Gauss-Kronrod quadrature scheme. We then compute the forward smile a t ,T with a 
simple root-finding algorithm. The Heston model parameters are given by p — —0.8, £ = 0.52, k = 1 and 
v = 6 = 0.07 unless otherwise stated in the figures. In Figures Q] and [2] we compare the true forward smile 
using Fourier inversion and the asymptotic in Theorem 14. II It is clear that the small-maturity asymptotic has 
very different features relative to "smoother" asymptotics derived in [3T]. This is due to the introduction of the 
forward time-scale and to the fact that the limiting mgf is not steep. Note also from Remark l4.2f iv) that the 
asymptotics in Theorem 14 . 1 1 can approach zero or infinity as the strike approaches at-the- money. This appears 
to be a fundamental feature of non-steep asymptotics; numerically this implies that the asymptotic may break 
down for strikes in a region around the at-the-money point. In Figure [3] we compare the true at-the-money 
forward volatility using Fourier inversion and the asymptotic in Lemma [4.3l Results are in line with expectations 
and the at-the-money asymptotic is more accurate than the out-of-the-money asymptotic. This is because the 
at-the-money forward volatility (unlike the out-of-the-money case) has a well defined limit as r tends to zero. 
In Figurc2]we use these results to gain intuition on how the Heston forward smile explodes for small maturities. 
In |31j the authors derived a diagonal small-maturity asymptotic expansion for the Heston forward smile valid 
for small forward start-dates and small maturities. In order for the small-maturity asymptotic in this paper to 
be useful, there needs to be a sufficient amount of variance of variance at the forward-start date. Practically 
this means that the asymptotic performs better as one increases the forward-start date. On the other hand 
the diagonal-small maturity asymptotic expansion is valid for small forward-start dates. In this sense these 
asymptotics complement each other. Figure [5] shows the consistency of these two results for small forward-start 
date and maturity. 




Lemma l4~3l and (|4.5[) conclude the proof. □ 
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FwdSmile Error 




Strike 



(a) Asymptotic vs Fourier inversion. (b) Errors. 

Figure 1. Here t = 1 and r = 1/24. In (a) circles, squares, diamonds and triangles repre- 
sent the zeroth, first, second and third-order asymptotics respectively and backwards triangles 
represent the true forward smile using Fourier inversion. In (b) we plot the errors. 



FwdSmile 

Error 




Snike 



(a) Asymptotic vs Fourier inversion. (b) Errors. 

Figure 2. Here t — 1 and r = 1/12. In (a) circles, squares, diamonds and triangles repre- 
sent the zeroth, first, second and third-order asymptotics respectively and backwards triangles 
represent the true forward smile using Fourier inversion. In (b) we plot the errors. 
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il 



ATMFwdVol 
0.28 r 




(a) ATM Asymptotic vs Fourier inversion. 



(b) Errors. 



Figure 3. Plot of the forward ATM volatility (r = 1/12) as a function of the forward-start 
date t. The Heston parameters arc p = —0.6, k = 1, £ = 0.4 and v = 6 = 0.07. In (a) circles, 
squares and diamonds are the zeroth-order, the first-order and the true forward ATM volatility. 



FwdSmik- 




t—'—t—i — i — i — i — 1 — i — i — i — i— i — i — i — i — i — 1 — i — i — l-1 — i — i — i — i—i — i — i — i — l_1 Strike 
L 0.8 0.9 1.0 1.1 1.2 1.3 

(a) t = 1, r = 1/100. 

Figure 4. Circles, squares, diamonds and tr 
third-order asymptotics respectively. 



FwdSmik' 




1 — i — , — , — u — i — i — i — L_. — i — i — i — — i — i — .—I — i — i — i — i — J — i — i — i — l_1 Strike 

0.8 0.9 1.0 1.1 1.2 1.3 

(b) t = 1, r = 1/1000. 
gles represent the zeroth, first, second and 



Fv, JSinilu 

0.50 r 

0.45 - 



0-94 0.96 0.98 1 .00 F02 1 .04 1 .06 

Figure 5. Here we compare the small- maturity third-order asymptotic (circles) to the diagonal 
small-maturity second-order asymptotic of [31] (squares) for t = 1/12 and r = 1/1000. 
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6. Proof of Theorems 13.11 and 14.11 

We split the prool of Theorems 13.11 and 14.11 into several parts, from Section 16.11 to Section 16.41 below. In 
Section 16.11 we develop the necessary tools to characterise the small- maturity Heston forward mgf domain and 
derive the Heston forward time-scale (Lemma 12. 3[) . In Section 16.21 we use the forward time-scale to define a 
time-dependent asymptotic measure-change and derive expansions for fundamental auxiliary functions needed 
in the analysis. In Section 16.31 we derive the asymptotics of the characteristic function of a re-scaled version 
of the forward price process (X^) under the asymptotic measure-change defined in Section 16.21 This section 
also uses Fourier transform methods to derive asymptotics of important expectations using this characteristic 
function expansion. Section 16.41 finally puts all the pieces together and proves Theorems 13.11 and 14.11 



6.1. Heston forward time-scale. A straightforward application of the tower property for expectations (see 
also [30) ) yields the Heston forward moment generating function: 



(6.1) A« (u, 1) := E (c^ t) )=A(u,t)+ vc~ Kt - ^ log (1 - 2f3 t B(u, r)) , for all u E V 

\ J 1 — 2p t B(u, t) £ J 

where 



t,TJ 



(6.2) 



. / k0 (, , .. „, / 1 — 7(w) exp (— d(u)T) 
A{u,t) := — [{n-p^u-d{u))T-2\og[ \ _ 

. k — p£u — d(u) 1 — exp (— d(u)T) 

B ( U ' T ) : = 72 1 r~\ — i At \ v 

5 1 — 7(u) exp (— d\u)T) 

1/2 K~ P iu-d{u) ? 



(6.3) d(u):=(( K -ptuf+u(l-u)e) , l(u):= K and [3 t := f (l - e""') . 

V / n — p^u + d[u) 4/t 

The first step is to characterise the forward time-scale in the Heston model. In order to achieve this we first 

need to understand the limiting behaviour of the re-scaled B function in (|6.2I) that plays a fundamental role 

in the analysis below. The following lemma shows that using h(r) = ^fr as a time scale produces the only 

non-trivial limit for the re-scaled B function. We then immediately prove Lemma 12.31 which characterises the 

forward time-scale in the Heston model. 

Lemma 6.1. Let h : R+ — > be a continuous function such that linii-x^o h(r) = and a £ M5_. The following 
limit then holds for B in (|6.2p for all u € R* : 

undefined, if t /h{r) /* oo, 

+00, if h(r) = or, 

+oo, if y/r/h^r) /■ oo and t /h(r) \ 0, 

0, if y/¥/h(r)\0, 

u 2 /(2a 2 ), if h^^ar 1 / 2 . 



limB(u/h(T),T) = < 



Proof. As t tends to zero we have the following asymptotic behaviours for the functions d and 7 defined in (|6.3D : 



(6.4) 



d (u/h(r)) = -L (K 2 h(r) 2 + uh(r) - 2np) - p^f 2 = j^d + d 1 + 0(h(r)), 
h{T) h(T) 

Kh(r) - p^u - ±ud - dxh{T) + O (h{T) 2 ) ih{r) ^ 

7 (u h{T)) = -—— ,/n tht\2\ = So .91 + O (h{r) ) , 

nh(T) — pt_u + iuo(o + dxhyr) + O (/i(r)^) u 



where we have set 

_-. , - t („\ a 1 ( 2k P - s § n ( u ) ip-psgn(u) (2k - p£) sgn(n) 

6.5 do := pCsgn(w), d x := — , g ~ — — 1—- gi := — —2- 

2p ip + psgn(w £p(p+ ipsgn(u)) 
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First let t/H{t) — > oo. Then exp (— d (u/h(r)) t) = exp (— iTp^\u\/h(r) + O (t)) , and so the limit is undefined 
(complex infinity). Next let t/H(t) = 1/a. Using (|6.4p we see that 

B(u/h(T),r) = (^P) ^ goe - imul/a + ° (1) = «/Mr) + O (1) , 

where C(u) := w (p£ c °t iu^p/2) — p£) , which is strictly positive for and ((0) = 0. It follows that the 

limit in this case is infinite. Next let t/H(t) — > 0. Here we can write 

B(u/h(r),r) = (-^±m + O (1)) ((-^r+O (k(r))) ( + O (r)) + O ^ T 



Zh{r) v 7 \\9o-l V W 7 V Kt) y ') \\h(r) 

fpU+±f)\u\\ ( 1 \ ±Tp£\u 



V £ J VflO-iy Kt) 



\2 



0{r/h{T)) 



2 

'r " 



= ~ 1 W) ) 



If y/r/h{T) tends to infinity, so does B(u/h(r), r). When ^pr jh(r) tends to zero then B(u/h(r), r) does as well. 
If • v /t//i(t) converges to a constant 1/a, then B(u/h(r),T) converges to u 2 /(2a 2 ), and the lemma follows. □ 

Proof of Lemma WM For any t > 0, the random variable Vt in (|2.2[) is distributed as /3 t times a non-central chi- 
square random variable with q = 4k8/£, 2 > degrees of freedom and non-centrality parameter A = ve~ Kt / fit > 0. 
It follows that the corresponding moment generating function is given by 

(6.7) A r( U ):=E(e^)=expf-^-)(l-2/3 tU )- 9/2 , for all u < -L. 

The re-normalised Hcston forward moment generating function is then computed as 
A«(u, h(r))/h(r) = E Je"^--**)/^)! = E [e ^(x t+T -x t )/h(r)^ 

= E ^(»/M^W+B(«/^WV t ) = e A(u/h(r),r) A V (fl ( tt /^( T ) ) T )) , 

which agrees with (|6.1[) when /i(t) = 1. This is only valid in some effective domain V t:T C R. The mgf for 
V* is well defined in V\ r := {u G R : B (u//i(t),t) < 1/(2^)}, and hence D t , T = D£ T n 2? r , where V T is the 
effective domain of the (spot) re-normalised Hcston mgf. Consider first T> T for small r. From [5J Proposition 
3.1] if £ 2 (u/h(T) - \)u/h{r) > (« - p^u/h(r)) 2 then the explosion time r*(u) := sup{i > : E(c" Xt ) < oo} of 
the Heston mgf is 

T* H (u/kM) = 2 :U 

Jf 2 (u/h(T\ - ThJhOA - (k- ntnlhirW 2 I 



7T ^{ P ^u/h(T)-K<0} 



Ve(u/h(r) - l)u/h(r) - (k - p£u/h(T)) 

Z ~ 7 \ \ o \ 

}■ 



arctan 



^/eWW) - Vu/h{r) - (k - ptu/h(T))* 



p^u/h(r) — k 

Recall the following Taylor series expansions, for x close to zero: 



arctan 



p^u/x — K 



\A 2 G _1 ) l~ ( k ~^) 2 ) = sgn(u) arctan (p) + ° (x) ' if p * °' 

arctan f — - A /^ 2 (- - l) - - K A =-I + 0(x), if p = 0. 
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As r tends to zero £ 2 (u//i(t) — l)u/h(r) > (k — p£u / 7i(r)) 2 is satisfied since p 2 < 1 and hence 

^7Tl {p=0} + \ (irl {pu < 0} + sgn(u) arctan (?) ) l Wo} + O (h{r))j , if u + 0, 

if u = 0. 



Therefore, for r small enough, we have t* h (u/h(r)) > t for all u £ R if r/h (r) tends to zero and (u/h{r)) > r 
for all u € (u_, if h (r) = ar, where 

2a /p\ 7ra 2a / f p\ \ 

u - ■■= ^ arctan ( - J !{p<o} - ^1{p=o} + ^- I arctan ^-J - tt J l {p>0} , 

2a f ( p\ \ ira 2a f p\ 

u + '■= I arctan I -J + tt J l {p<0} + ^-l{ P =o} + arctan ^- J l {p>0} . 

If t//i(t) tends to infinity, then (u/h(r)) < r for all !ieK*. We are also required to find T>Y T for small r. 
Using Lemma 16.11 we see that if h(r) = ar 1 / 2 then lim rV) l>tf T = {" G K : \u\ < a/y/JT t }. If T^ 2 /h(r) tends 
to infinity then lim T ^o T>Y T = {0} and if it tends to zero, then lim^o T^Yr = The limiting domains in the 
lemma follow after taking the appropriate intersections. Next we move on to the limits. We only consider the 
cases where h(r) = ar 1 / 2 and where t 1 / 2 /K{t) tends to zero since these are the only cases for which the forward 
moment generating function is defined. Using (|6.6|) we see as r tends to zero 

!RSi s , h MR , ,, n ,v B(u/h( T ),T)ve- Kt $ 0(1), if = arVa, 

6.8 log 1 - 2a t S (u/h(T), r) = - . — — — - — - = < 

1 ' ^ ' K> " l-2p t B{u/h{T),r) \ 0(T/h(r)), if y/¥/h(r) \ 0. 

The lemma follows from this and the fact that the function A in (|6.2[) satisfies A(u/h(T),r) = O ((t//i(t)) 2 ) . □ 

6.2. Asymptotic time-dependent measure-change. In this section we define the fundamental asymptotic 
time-dependent measure-change in (|6.15j) and derive expansions for critical functions related to this measure- 
change. In order to proceed with this program we first need to prove some technical lemmas. We use our 
forward time-scale and define the following rescaled quantities: 

(6.9) AM(u):=A?\u,V^), A{u):=A(u/y/r,r), B(u) := B(u/V¥,t), 

with At\ A and B defined in (|6.1j) and (|6.2j) respectively. The following lemma gives the asymptotics of the 
re-scaled quantities A, B as r tends to zero: 

Lemma 6.2. The following expansions hold for all u £ T>\ as r tends to zero (B\ was defined in (|3.2p ): 

_ _ _ ii 2 kQi~ 

(6.10) B(u) = y + B 1 (u)^ + 0(t) 1 A(u) = ^ r + (D(T 3 / 2 ). 

Proof. From the definition of A in (|6.2[) and the asymptotics in (|6.4[) with h(r) = ypr we obtain 

m := A r) - * ((. - f - , (i)) T - 21 „ g (l^tefc^)) 



k6 f f p^u ±udo 
-2 log 



-dl + 0(V7) T 



' 1 - (. 9 o - i-s/r^/u + 0(t)) exp (-iud V^ - diT + 0(r 3 / 2 )) ' 



l-(.g -iV^.9i/« + O(r)) 
u 2 6»kt/4 + C(t 3 / 2 ). 
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Substituting the asymptotics for d and 7 in (|6.4j) we further obtain 

1 - exp (-d(u/y/T)T) 1 - exp (-±ud ^/T - dir + C(t 3 / 2 )) 



1 - 7 ( U /Vr)exp (-d(«/Vr)T) 1 - (,g - iVrfli/u + 0{r)) exp (-i<zd ^ - diT + 0{t*/ 2 )) ' 
and therefore using the definition of B in (|6.2j) we obtain 



.12) B(u):=B \-=,r = 



k - pju/y/r - d{u/y/r) 1 - exp (— (u/yff) r) 

£ 2 1 - 7 (w/V? ) exp (-d (u/y/r) r) 



+ Si(u)VF+ 0(r) = — + B 1 (u)VT + 0(t). 
S 1 - .90 2 

□ 

It is still not clear what benefit the forward time-scale has given us since the limiting mgf is still degenerate. 
Firstly, even though the limiting mgf is zero on a bounded interval, the re-scaled forward mgf for fixed r > is 
still steep on the domain of definition which implies the existence of a unique solution u*(k) to the equation 

(6.13) d u A®(u*{k))=k. 

Further as r tends to zero, u*(k) converges to l/y/Pt when k > and to —l/yffit when k < (see Lemma 
below). The key observation is that the forward time-scale ensures finite boundary points for the effective 
domain, which in turn implies finite limits for u*(k). This is critical to the asymptotic analysis that follows and 
it will become clear that if any other time-scale were to be used the analysis would break down. The following 
lemma shows that our definition (|6.13[) of u*(k) is exactly what we need to conduct an asymptotic analysis in 
this degenerate case. 

Lemma 6.3. For any k G K, the equation (|6.13p admits a unique solution u* T (k) which converges as r tends to 
zero to 1/ ' \/pt (—\J\fpt) when k > (k < 0) and to zero when k = 0. Further u*(k) G T>\ for r small enough. 

Proof. From Lemmas 12.31 and 16.21 for u G T>\ := (— l/y 7 /^", l/y 7 /^) and r small enough, we have 

(6.14) d u A^{u)=A (u)V : F + <D(r), 
where we set 

4ku (9e 2Kt (4k - £_ 2 u 2 ) + 2e Kt (0£ 2 u 2 - 2k6 + 2nv) - 9£ 2 u 2 ) 



(c Kt (4k - £ 2 u 2 ) + £, 2 u 2 y 
The following statements are straightforward (Figure [5] provides a visual help) : 

(i) Ao is a strictly increasing map from l/y/JJ^) to R; 

(ii) Ao(0) = 0. 

Consider the case k > (the case k < follows by symmetry). For sufficiently small r, u*(k) G T>\, since the 
expansion in (|6.14j) holds on 2?^, and (i) implies that a unique solution to (|6.13j) exists. For small enough r, 
u*(0) = from (ii), i.e. the origin is a fixed point. The function Ao is strictly positive on M5_, and hence for 
small enough t, u*(k) is strictly increasing and bounded above by l/y/ftt, and therefore converges to a limit 
L G [0, l/y/pt\. If L G [0, then the continuity of Ao implies that lim r x^o Ao(u* (k))y/r = 7^ k, which is 
a contradiction. It follows that L = l/y/Jh and the lemma follows. □ 



For sufficiently small r we introduce a time-dependent change of measure by 
3.15) ^ :=cxp(<(fc)xW/V7- AW«(fc))/V7) . 
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FIGURE 6. Plot of u H> Ao(u)y / r for different values of r. Circles, squares and diamonds 
represent r = 1/2, 1/12, 1/50. The forward-start date is t = 1 and the Heston model parameters 
are v = 9 = 0.07, £ = 0.4, p = -0.6. The limiting domain is VV^t) ~ (-6.29, 6.29). 

By Lemma [6~3l u*(k) £ T>\ for r small enough and so jAr^u*)! is finite since T>\ = liuir^olu £ K : |Ar \u)\ < 
oo}. Also dQfe T /dP is almost surely strictly positive and by definition E[dQ^ T /dP] = 1. Therefore (|6.15[) is a 
valid measure change for all k £ M* and sufficiently small r. Equation (|6.13[) can be written explicitly as 

( 6 - 16 ) ^^[e^ t A\ul)(A-2B{ul)p t )\B\ul)( K AKep tC ^ = (l - 2B{u* T )p t )\ 

with A and B defined in (|6.9[) . We now use this to derive an asymptotic expansion for u* as r tends to zero. 

Lemma 6.4. The expansion u*(k) = a a (k) + ai(fc)r 1 / 4 + a 2 (fc)r 1 / 2 + a 3 (fc)r 3 / 4 + O(t) holds for all k £ R* as 
t tends to zero, with ao, a±, a 2 and 03 defined in (13. ip . 



Proof. Existence and uniqueness was proved in Lcmma l6.3l and so we assume the result as an ansatz. Consider 
k > 0. From Lemma 16.31 it is clear that ao(fc) = l/y^t- The ansatz and Lemma 16.21 then imply the following 
asymptotics as r tends to zero (we drop here the fc-dependence): 

B(u*) =^jf t + aoaiT 1/4 + rr 1 / 2 + (ena 2 + a a 3 + aiB[(a )) r 3 / 4 + O(r), 

(6.17) B'[u*) =a + air 1 / 4 + (a 2 + ^ (ao))^/ 2 + G(r 3 / 4 ), 

=i^a r + O(r 5 / 4 ), 

where r = r(k) := a a 2 + -Bi(a ) + a 2 /2 = af/2 — nd/ {\k\£? \fj3t) is defined in (|3.3[) . We substitute these asymp- 
totics into (|6.16j) and solve at each order. At the r 1 / 4 order we have two solutions, a\ (k) = ±^/vc~ Kt / 2 / (2\/kf3^ 4 ) 
and we choose the negative root so that u* £ 2?^ for r small enough. In a straightforward, yet tedious, manner 
we continue the procedure and iteratively solve at each order (the next two equations are linear in a 2 and 03) 
to derive the asymptotic expansions in the lemma. An analogous treatment holds in the case k < 0, and the 
lemma follows. □ 

In the forthcoming analysis we will be interested in the asymptotics of 

(6.18) e T (k) := (l - 2B«(Af))ft) r" 1 / 4 , 

as t tends to zero. Since (1 — 2B(u* (k)){3 t ) converges to zero, it is not immediately clear that e T has a well-defined 
limit. But we can use the asymptotics in f|6 . 1 T[) to deduce the following lemma: 
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Lemma 6.5. The expansion e T (k) = eo(fc) + ei(fc)r 1 / 4 + e2(fc)r 1 / 2 + C(r 3 / 4 ) holds for all k £ M* as r tends 
to zero, where eo, ei and e2 are defined in 



Proof. We substitute the asymptotics for B(u*) in (|6.17[) into the definition of e T in t)6. 18[) and the lemma 
follows after simplification. □ 

6.3. Characteristic function asymptotics. We now define the random variable Z Tt k ■= (x^ — kj /r 1 / 8 

and the characteristic function : M. — >• C of Z T ^u in the Qfc jT -mcasurc in (|6.15[) as 



(6.19) 

Define now the functions c\>\ , i 



4>i(k,u) := iu 7p (k) + 



4a o (fc)0 K /3 t 



$ T ,k(u) :=E Qfc - (e iuZ - fc ). 
^Cby 

+ iu 3 V'i(fc), and ^ 2 (M) := u 2 </)%(k) + u A (f> b 2 {k) + u 6 <p c 2 (k), 



eo(k)e 

with -00; 0i defined in (|3.5[) , ao, eo in (|3.1[) , p.4[) . and ^f; 02 an< ^ 02 m (|3.6I) . The following lemma provides 
the asymptotics of & T ,k'- 

Lemma 6.6. The following expansion holds for all k £ R* as r tends to zero (with a given in (|3.3p ): 
$r,fc(«) - c"^" 2 / 2 (l + ^(jfe, ^t 1 / 8 + 2 (fc, u)r 1 /4 + o (r 3 / 8 )) . 

Remark 6.7. For any tgl', Levy's Convergence Theorem [57J Page 185, Theorem 18.1] implies that Z T ^ 
converges weakly to a normal random variable with zero mean and variance cr 2 (k) as r tends to zero. 

Proof. From the change of measure (|6.15[) and the forward moment generating function given in (|6.1[) we 
compute (we drop the fc-dcpendcncc throughout) for small r: 



log$ fc ,r(w) = log I 



= logE p 



exp 



i*{k)xi t] a^co' 



exp 



luXi^ luk 



iA?)«)-^ + lo g E' 



exp 



(±ut 3/s + u*\ 



(6-20) = - ™ + -1 (A« (i^r 3 / 8 + <) - A« «)) . 

Using the asymptotics in (|6.17[) we have that as r tends to zero (we drop the fc-dependence) 
(6.21) B (< + i«r 3 / 8 ) = ^ + aoaxr 1 / 4 + ia r 3 / 8 u + rr 1 / 2 + G(r 5 / 8 ), 



B (<) = ^ + aoair 1 / 4 + rr 1 / 2 + C(r 3 / 4 ), 
B (< + iur 3 / 8 ) - B «) = ±a Q uT 3 ^ + ia lU T 5 / 8 - \u 2 r 3 ^ + 0(r 7 / 8 ), 
where r = r(k) := a a 2 + B 1 (ao) + a 2 /2 = a 2 /2 — k6 / (\k\^ 2 y/]3t) is defined in (|3.3[) . Similarly for small r, 



3.22) 



A (< + i^r 3 / 8 ) = ^a 2 r + ^a oai r 5 / 4 + ^a,ur^ + 0{r 3 ' 2 ) 



l( <) = ^ a 2 r+ ^ aoaiT 5/4 + 0(r 3/ 2)) 



A (< + inr 3 / 8 ) - 1«) = ^-aour 11 / 8 + 0(t 3 / 2 ). 
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We now use e T denned in (|6.18p to re-write the term 

B (u* + iur 3 / 8 ) ve~ Kt vQ- Kt T- x ^B (u% + iur 3 / 8 ) 



1 - 2p t B (u* + iur 3 / 8 ) 6t _ 2(3 tT - 1 / 4 (S ( u* + iur 3 / 8 ) — B (u*) 

and then use the asymptotics in (|6.2ip and Lemma 16.51 to find that for small r 
ve - K t T -i/4g (u* + iur 3 / 8 ) 



(6.23) 



e r - 2/3 t r-V4 (b (u* + iur 3 / 8 ) - B (u* T ) 

ue -«V-V4 (al/2 + aoair 1 / 4 + ia r 3 / 8 u + rr 1 / 2 + 0(r 5 / 8 )) 
~ e + dri/4 + e2r i/2 + 0( r 3/4) _ 2f3 tT - 1 / 4 (ia ur 3 / 8 + iaiur 5 / 8 - ±u 2 r 3 / 4 + C(r 7 /8)) 

W e- K ^a| 1/4 ue~ Bt ia 3 u/3 t i/s , -«t /" £o£i _ ogei \ 



2e eg V e 2eg / 

2„.2 

+ (iu^o + iu 3 ^'i)r 1/8 + (V>4 + ^2U 2 + V^ 3 w 4 )t 1/4 + G(r 3 / 8 ), 



2 

with cr and tpo , . . . , "04 defined in (|3.3p and (|3.5p . From the definition of do, eo and /?t we note the simplification 

i^e~ Kt a^(fc)u/3 t _ iufc 
1 J egCJfcJrVs T i/8- 

Similarly we find that as r tends to zero 

B «) we- K * vc-^t- 1 / 4 ^ «) _ we-^r" 1 / 4 (a 2 /2 + ao^r 1 /* + rr i/2 + 0( T 3/4)) 



(6.25) 



l-2/3 t B«) e r e + eir 1 /^ + ear 1 ^ + 0( r 3/4) 



-,2,,„-Kt 



2e V e o 2 eo / 

Again we use e T defined in (|6.18p to re- write the term 



exp 



2k6 , / 1 - 2^5 (u*) 
log 1 



u* + iur 3 / 8 ) 



= 1- 



2/3 t (B « + iur 3 / 8 ) -B(u* T ) 
fiTrV 4 



and then use the asymptotics in (|6.2ip and Lemma 16.51 to find that for small r 

\ _ ^(^< + ^p-B {K) ) ^ ~ ¥ = ^ + 2_^U t1/8 4a^ 1/4 + 0(t3/8) ^ ^ 

6.26) =1 + ^-^r 1/8 — v 9 IJ-T 1 ' 4 + C(r 3/8 . 

eoC £ eg 

Using (|6~2^1) with definition (EHI) and (|63|) . property (|(T2"3)) and the asymptotics in (|rT2"2"]) , (|6T2"5j) 
and (|6.26p we finally calculate the characteristic function for small r as 

- + (iu^o + iu^r 1 / 8 + (fou 2 + ^u 4 )t 1 I 4 + G(r 3 / 8 )^) 
' 4i K 6a p t u 1/8 4^a 2 /3 2 u 2 (2^ + £ 2 ) 1/4 3/ \ 

e e 2 ^ +CAr >r 



^k,r(u) = exp - 
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with tpo , . . . , tp3 defined in Q3.3P , (|3.5[) , and so the lemma follows from the following decomposition 
2 - 2 \ ' 4K0ao/3t 



r(u) = cxp ^-J |l + i + 



e £ 2 



.1/8 



V>g 4K<9a 2 ,/3 2 u 2 (2k(9 + £ 2 ) 4K6»a /3 t 



e £ 2 



-0o u + "03 - 0o"0i - 0: 



e e 



4 wVi 



■*+0(r*)}. 



□ 



The following technical lemma will be needed in Section 16.41 where it will be used to give the leading order 
exponential decay of out-of-the-money forward-start options as r tends to zero. 

Lemma 6.8. The following expansion holds for all k € R* as r tends to zero: 

e -K/v^+A?'«)/V? = c -^(k)/^+c o (k)/r 1 ' 1 +c 1 (k) T -K0/(2e) C2 (J^ + Zl (fc) T V4 + 0(t 1/2 ) \ , 

where cq, c\ and c-i are defined in (|3.8p . A* is characterised explicitly in Lemma\2.J\ and Z\ is given in (|3 . T[) . 



Proof. We use the asymptotics in Lemma 16.41 and the characterisation of A* in Lemma 12.41 to write for small r 
(we drop the fc-dependence) 



(6.27) 



exp 



(-ku*/y/r) = cxp (-a k/Vr - aifc/r 1/4 - a 2 fc) (l - a 3 kr 1/4 + C(r 1/2 )) 

(-A*(jfe)/VT - aifc/^ 1/4 - a 2 fc) (l - agfcr 1 / 4 + ©(r 1 / 2 )) 



cxp 



Using the Hcston forward mgf definition in (|6.1[) and (|6.9[) we can write 

B{u* T )ve~ Kt 



(6.28) 



exp 



(AW«)/V^) =exp U«) + 



^log(l-2A-B«))l ■ 



1 - 2[3 t B(u* T ) e 

Using the definition of e T in f|6. 18[) and the asymptotics in Lemma 16.51 we find that for small r 
(6.29) exp (~iog(l - 2&B«))) = T -"*/W)e;™/? = T -^/( 2 ? 2 ) (e + + ©(r 1 / 2 )) 



_ - K e/(2e) p -^s/e 



1 _Mei T L,* + 0{T L,2 ) 
£ e 



The the lemma follows after using (j6~2"7j) and (j6~2"8"]l . the asymptotics in (|6~29"]) . (|6T2"5j) and ((6~2"2"j) and the 

simplification c (A;) = ve~ Kt /(2e (k)[3 t ) — a\(k)k = 2\a\(k)k\. □ 

We now use the characteristic function expansion in Lemma 16.61 and Fourier transform methods to derive 
the asymptotics for the expectation (under the measure (|6.15l) ) of the modified payoff on the re-scaled forward 
price process. This lemma will be critical for the analysis in Section [ 



Hfc<o} 



Lemma 6.9. The following expansion holds for all k € R* as r tends to zero: 



l {fe>0} 



(l+ Pl (ky/ 4 + o(r^)), 



a(k)V2TT 

where a is defined in (|3.3|) . p\ in (|3.7p and j3 t in (|6.3p . 
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Proof. We first consider k > and drop the /c-dependence for the functions below. We denote the Fourier 
transform T by (Ff)(u) := e lux f(x)dx, for all / € L 2 , u £ R. The Fourier transform of the payoff 



-<Zr, k /T 3/a ( pZ^kT 1 '* _ 1 \ ig given by 



(iU-U*/T 3 / S +T 1 / 8 ) 
1 



eZ (i«-u;/r 3 / 8 ) 

(iu — it*/r 3 / 8 ) 



1 



" (ill - u*/r 3 / 8 ) (iu - u*/r 3 / 8 + r 1 ^) 

r 7/8 

(u* - ir 3 / 8 u) (u* — y/r- ir 3 / 8 w) ' 

if u* > max(T 1 / 2 ,0) = t 1 / 2 , which holds for r small enough since u* converges to a > by Lemma \6. 41 Due 
to Remark 16.71 Z T converges weakly to a Gaussian random variable and since the Gaussian density and the 
modified payoff are in L 2 we can use Parseval's Theorem [25) Page 48, Theorem 13E] for small enough t to write 

r 7 / 8 $ r , fe (u) 



(6.30) 

where we have used that 



+ iT 3 / 8 it) (u* — s/t + iT 3 / 8 w) 



du, 



r7/8 



r 7/8 



(u* — ir 3 / 8 u) (u* — s/t — ir 3 / 8 u) (u* + ir 3 / 8 u) (u* — ^fr + ir 3 / 8 u) ' 

with a denoting the complex conjugate for a G C. Using the asymptotics of u* given in Lemma 16.41 we can 
Taylor expand for small t to find that 

T 7/8 T 7/8 



,7/8 



3.31) 



i _ + L 

an V 



.3/8 



' (u* + ir 3 / 8 u) (u* - V? + ir 3 / 8 u) a 2 + 2a a 1 T 1 / 4 + C(t 3 / 8 ) a 2 \ a 

Finally combining (|6.3ip and the asymptotics of the characteristic function derived in Lemma 16.61 with (|6.30p 
we find that for small r 

l_ f°° tV^t^u) 

2tt 



VZ + 1T 3 / 8 U I I UZ. — \/T + 1T°' °U 



3/8, 



-du 



r 7/8 



e 2 



(l + 0i (u, fc)r 1 / 8 + ^ 2 (u, k) - ^) + o ( r 3 / 8 ) ^ du 



T 7/8 /-00 2 2 

e"^ - ( 1 + ( u 2 d>% + u A <b\ + u 6 - 



a^n J _ 00 



^l)r^ + 0(r"/ s ))du, 



do 



where in the last line we have used that e -£ ~r !- (/>i(u, k)du = 0, since <j)\ is an odd power of u. The result 
then follows by using simple moment formulae of the normal distribution. Fix now k < 0. The Fourier transform 
of the payoff q- u *t z ^^I t ^ i& [1 — c ZT kTl/8 j is given by 



"dz 



(iu - w*/r 3 / 8 ) 
1 

(iu — u*/r 3 / 8 ) (iu — u*/t 3 / 8 + r 1 / 8 ) 

^7/8 



(iu-^/r^ + T 1 / 8 ) 
1 



i* — iT 3 / 8 w) (u* — s/t — It 



3/8, 
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if it* < min(T 1 / 2 , 0) = 0, which holds for r small enough since u* converges to ao < by Lemma T6. 41 The rest 
of the proof is analogous to k < above and we omit it for brevity. □ 

Remark 6.10. We have chosen to specify the remainder in the form o(l/r 1 / 4 ) instead of C(l/r 3 / 8 ) since it 
can actually be shown that the term 0(l/r 3 / 8 ) is zero by extending the results in Lemma [6.61 and the next 
non-trivial term is ©(I/t 1 / 2 ). For brevity we omit this analysis. 

6.4. Option price and forward smile asymptotics. In this section we finally put all the pieces together 
from Sections 16.11 - 16.31 and prove Theorems 13.11 and 14.11 

Proof of Theorem \3.1\ We use the time-dependent change of measure defined in (|6.15|) to write forward-start 
call option prices for all k > as 



E 



3 A««)/v^ F Q fe , 



k-4"(.;i 

y/7 



e ^ 



(e*—e*)" 
^(e^ V8 -l)' 



with Z T k defined on page 1171 A similar result holds for forward-start put option prices for all k < 0. The 
theorem then follows by applying Lemma 16.81 and Lemma 16.91 and using put-call parity since in the Hcston 
model (e Xt )t>o is a true martingale j4] Proposition 2.5]. □ 



Proof of Theorem \4-l\ The general machinery to translate option price asymptotics into implied volatility 
asymptotics has been fully developed by Gao and Lee [22]. We simply outline the main steps here. Assume the 
following ansatz for the forward implied volatility as r tends to zero: 

°IM = + ^ + V2{k , t ) + v s (k, tyi* + o(r^). 

Substituting this ansatz into the BSM asymptotics in Lemma 13.41 we then obtain 



exp 



k 2 k 2 v\ 




1 



k 2 (vf - 2u wiw 2 + v 2 v 3 ) 



'2-Kk 2 

Equating orders with Theorem 13.11 we solve for vq and v\ , but we can only solve for higher order terms if 

T 3/4 = r (7/8-e«/(2 f 2 )) Qr 4k0 = n 
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